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A three-dimensional ±J XY spin-glass model is investigated by a nonequilibrium relaxation 
method. We have introduced a new criterion for the hnite-time scaling analysis. A transition 
temperature is obtained by a crossing point of obtained data. The scaling analysis on the relax- 
ation functions of the spin-glass susceptibility and the chiral-glass susceptibility shows that both 
transitions occur simultaneously. The result is checked by relaxation functions of the Binder pa- 
rameters and the glass correlation lengths of the spin and the chirality. Every result is consistent if 
we consider that the transition is driven by the spin degrees of freedom. 
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I. INTRODUCTION 

Random and/or frustrated systems exhibit a variety of 
exotic phenomena. Common knowledge of uniform sys- 
tems may not be applied to these systems. Randomness 
sometimes brings order to disorder. Exotic states may 
appear in the ground state. New-type phase transitions 
occur in these systems. They are referred to as complex 
systems and have been attracting wide interest. 

Spin glasses^ are prototype of complex systems. Fer- 
romagnetic interaction and antiferromagnctic interaction 
are randomly distributed in the materials. There is frus- 
tration in local interaction bonds. Magnetic spins cannot 
always minimize the interaction energy. They freeze at 
low temperatures while pointing random directions. This 
freezing is called the spin-glass transition. 

Theoretical investigations on spin glasses are mainly 
made into the Edwards- Anderson models. Spins are lo- 
cated on a regular lattice but the interaction bonds be- 
tween spins are randomly distributed. Many theoreti- 
cal methods have been developed. By progress of com- 
putational facilities numerical simulations are recently 
performed. However, simulations in spin glasses suffer 
from serious slow dynamics. It takes a long time to re- 
alize equilibrium states at low temperatures. A numer- 
ical method overcoming this difficulty can be success- 
fully applied to other complex systems. One example is 
an application of the temperature-exchange Monte Carlo 
method 3 - to protein folding problems^. 

Not many problems have been made clear due to the 
difficulty of simulations in spin glasses. One unsolved 
problem is whether or not a theoretical model with 
short-range interaction bonds can explain the real spin- 
glass transition. There is a general consensus that the 
spin-glass transition occurs in the Ising model in three 
dimensions The estimated critical exponents agree 
with the corresponding experimental results^. Spins of 
many spin-glass materials are of Heisenberg type. The 
Heisenberg spin-glass model in three dimensions should 
exhibit the spin-glass transition. However, numerical 
investigations2*i2*ii concluded that there is no spin-glass 
transition in this model. There has been discrepancy be- 
tween theory and experiment in this model. 

There are three major theories explaining this dis- 



crepancy. One is the anisotropy theory. There are fi- 
nite magnetic anisotropics in exchange interactions in 
real materials. Even though there is no spin-glass order 
in the isotropic Heisenberg model, the spin-glass transi- 
tions in real materials are explained by finite anisotropy. 
Matsubara et al. showed that the anisotropy reinforces 
the spin-glass orderii 2 . The second theory is the chiral- 
ity mechanism introduced by Kawamurapi 3 ^ It is ar- 
gued that there occurs the chiral-glass transition with- 
out the spin-glass transition. Finite anisotropy in real 
materials mixes the chiral degrees of freedom and the 
spin degrees of freedom. This mixture links the chirality 
and the spin. The spin-glass transitions in real materials 
are driven by this chiral-glass transition. This scenario 
is based on that there is no spin-glass transition in the 
isotropic model. Since the chirality is defined by spin 
variables, the chiral-glass transition trivially occurs if the 
spin-glass transition occurs. Therefore, it is crucial to 
check whether both transitions occurs at the same tem- 
perature or not. The third theory argues a possibility 
of this simultaneous transition. Matsubara et alMti^iL 
showed several lines of evidence that the spin-glass transi- 
tion occurs at a hnite temperature. Nakamura et 
investigated the model by a nonequilibrium relaxation 
mcthod22i2ii22i2i24 and obtained the consistent results. 
Here, we note that nonequilibrium relaxation of the spin- 
glass susceptibility 2 ^ and that of the Binder parameter 
are known to exhibit the algebraic divergence at the spin- 
glass transition temperature in the Ising model. Lee and 
Young2£ also concluded a single chiral-glass and spin- 
glass transition in the model with Gaussian random bond 
distributions. These results suggest that the previous 
investigations which concluded no spin-glass transition 
have some technical problems. 

The situation is same in the XY spin-glass model in 
three dimensions. Spins are restricted to point in the XY 
plane in this model. A domain- wall energy analysis^i and 
a Monte Carlo simulation analysis^ concluded no spin- 
glass transition. The chirality mechanism was originally 
proposed in this modeli22i22*2I However, a possibility of 
the finite spin-glass transition has been suggested by re- 
cent investigationa 27 i 32 i 33 ' 34 i 35 . 

In this paper we focus on the problem in the XY spin- 
glass model in three dimensions: whether or not the spin- 
glass transition and the chiral-glass transition occur at 
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the same temperature. It is also made clear why there has 
been the discrepancy even though investigating the same 
theoretical model with similar methods. Wc consider the 
following three points as origins of the discrepancy. One 
is ambiguity in the scaling analysis. The obtained results 
are sometimes dependent on a temperature range of data 
used in the analysis. We introduce a new scaling criterion 
to eliminate this ambiguity. The second point is strong 
effects from system sizes and boundary conditions. We 
carefully observe how these effects appear in the relax- 
ation data. It is found that the effects are very strong and 
strange. The third point is a use of the spin-glass sus- 
ceptibility and the Binder parameter in the equilibrium 
simulations on finite lattices. Both quantities are under 
the strong influence of the size effects. These points are 
discussed in this paper. 

The present paper is organized as follows. Section [n] 
describes the present model and our method. A new 
scaling procedure is explained with an example of the 
ferromagnetic Ising model in two dimensions. Section ITTT1 
shows our results. A finite-time-scaling analysis on the 
glass susceptibility is performed. The results arc con- 
firmed by relaxation functions of the Binder parameter 
and the glass correlation lengths. Section llVl is devoted 
to summary and discussions. 



is defined by the following expression. 



Xsg 



i 

N 



h3 



Si • S j 



2N 



m(m — I) 



E (O 2 

A>B,fi,v 



(2) 

The subscripts /j, and v denote two components, x and y, 
of spin variables S. The thermal average is denoted by 
(■■•), and the random-bond configurational average is de- 
noted by [• • • ] . The thermal average is replaced by an av- 
erage over independent real replicas: (• • ■ ) — > (1/to) Y^a- 
Replica number is denoted by to. Superscripts A and B 
are the replica indices. The number to determines reso- 
lution of the thermal average. It should be large in order 
to improve the accuracy of average. We prepare 32 repli- 
cas (to = 32) with different initial spin configurations in 
this paper. Each replica is updated in parallel with a 
different random number sequence. An overlap between 
two replicas is defined by 
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The chiral-glass susceptibility^ is similarly defined by 
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II. MODEL AND METHOD 
A. Model and physical quantities 

A model we consider in this paper is the ±J random 
bond XY model in a simple cubic lattice. The Hamilto- 
nian is given by 
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where the sum runs over all the nearest-neighbor spin 
pairs (ij). Each spin, S, is defined by angle, 0j, in the XY 
plane: Si = (cos0j,sin#j). The angle 6i takes continuous 
values in [0, 2ir), but we divide it to 1024 discrete states 
in this paper. Effects of the discreteness are checked to 
be negligible by comparing data of 1024-state simulations 
with those of 2048-state simulations. The interactions Jy 
take two values of + J and — J with the same probability. 
The temperature T is scaled by J . Linear lattice sizes 
are denoted by L. Total number of spins is N = L x 
L X (L + 1), and skewed periodic boundary conditions 
are imposed. Spins are updated by the single-spin-flip 
algorithm using the conventional Metropolis probability. 

Physical quantities observed in our simulations are the 
spin-glass susceptibility Xsg ; the chiral-glass susceptibil- 
ity Xcg; the Binder parameters g sg and g cg , the glass cor- 
relation length £ sg and £ cg . The spin- glass susceptibility 
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A local chirality variable, . IS defined at each square 
plaquette a that consists of the nearest-neighbor bonds 
(Jk)- 

The Binder parameters for the spin glass and for the 
chiral glass are calculated through the following expres- 
sions. Here, we have also replaced the thermal averages 
by the average over real replicas. 



= 3-2 



E 



m(m — l) 

A>B 



\ 



E 



2 ( n AB)2 

m(m—l) £-i Wfi" I 
y ' A>B 



3c g 



i(m — l) 



IT E (<z2 



AB\i 



A>B 



\ 



i(m — l) 



it E (q AB Y 



A>B 



(8) 



Since the Binder parameters lose size dependences at the 
transition temperature, the noncquilibrium relaxation 
functions exhibit algebraic divergence with an exponent 
d/zm^ From the exponent we can obtain a value of the 
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dynamic exponent z. It also serves as another check of 
the spin-glass transition. If the susceptibility and the 
Binder parameter exhibit algebraic divergence, it is very 
certain that the phase transition occurs. Kawamura and 
Li2l argued that the Binder parameter of the chiral glass 
shows a negative dip at the transition temperature. They 
observed crossing behavior of the Binder parameter on 
the negative side. We observe noncquilibrium relaxation 
of the Binder parameter in order to check this negative 
crossing: it should diverge algebraically on the negative 
side. 

A spin-glass correlation function is defined by the fol- 
lowing expression. 
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A length between two spins is denoted by I. Here, 
q^ B (i) = s\ A ^ s\ B J is a spin overlap on the i site. A 
chiral-glass correlation function is defined by the follow- 
ing. 



/cg(0 



- V V 

3m(m - 1) ^ ^ 

v ' A>B a 



q^ B (a)q^(a + l) 



(10) 



k1 A 'k1 B ' is a local overlap of the chiral- 
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ity. We only consider a distance / between two parallel 
square plaquettes. We estimate a spin-glass correlation 
length £ sg and a chiral-glass correlation length £ cg fitting 
the correlation functions by the following expression. 
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algebraically as £ S g/cg oc t 1 /* . We try to find such a tem- 
perature first observing raw relaxation functions. If a 
temperature is higher than the transition temperature, 
the susceptibility converges to a finite value even in the 
infinite-size system. The lowest temperature at which 
the relaxation data exhibit the converging behavior is 
the upper bound for the transition temperature. 

The second method to estimate the transition temper- 
ature is the finite-time scaling analysis^. This is a di- 
rect interpretation of the conventional finite-size scaling 
analysis by considering the dynamic scaling hypothesis: 
t ~ £ z . A term "size" is replaced by a term "time" 
by this relation. We collect relaxation functions of the 
spin-glass susceptibility at various temperatures. They 
should fall onto a single universal curve if we properly 
choose the transition temperature, T sg , and critical ex- 
ponents, 7 and zv. The following equation is the scaling 
formula. 

Xs e (t)t- 7/zV = Xs S (t/r(T)), t(T) « (T - T sg )— (12) 

Here, x denotes a universal scaling function. Monte Carlo 
step t is scaled by correlation time t(T) which diverges 
algebraically at the transition temperature. We plot 
data of the spin-glass susceptibility multiplied by t^ 1 ! 2 - 11 
against t/(T — T sg )~ 2iy . An estimation of the chiral-glass 
transition temperature is performed by the same proce- 
dure. 

At the estimated transition temperature we observe re- 
laxation functions of the Binder parameter and the corre- 
lation length. They should exhibit algebraic divergence. 
The observation supports the occurrence of the transi- 
tion. The dynamic exponent z is estimated by these re- 
laxation functions. The exponent of the Binder parame- 
ter is d/z and that of the correlation length is 1/z. We 
can check consistency of the estimates comparing these 
two values. 



Exponent d is a dimension of the lattice [d = 3), 7 is the 
glass susceptibility exponent and v is the glass correlation 
length exponent. 



C. New criterion for scaling analysis 



B. Nonequilibrium relaxation method 

We start simulations from random spin configurations. 
This is a paramagnetic state at T = 00. The tempera- 
ture is quenched to a finite value from the first Monte 
Carlo step. At each step we observe physical quantities 
and obtain the relaxation functions. Changing an initial 
spin state, a random bond configuration, and a random 
number sequence we start another simulation and ob- 
tain another set of the relaxation functions. We take 
averages of relaxation functions over these independent 
Monte Carlo simulations. The relaxation functions of 
the susceptibility and the Binder parameter multiplied 
by N exhibit algebraic divergence at the transition tem- 
perature: Xsg/cg t~</ zy and g S g/ C g x N oc t d l z . A re- 
laxation function of the correlation length also diverges 



Finite-size/time-scaling analysis is a powerful tool to 
investigate critical phenomena. The method has been ap- 
plied to various systems successfully. However, it some- 
times has ambiguity obtaining correct values of the tran- 
sition temperature and the critical exponents. One exam- 
ple is a spin-glass transition problem in the ± J Hciscn- 
berg model in three dimensions. Olive et aliS> performed 
a finitc-size-scaling analysis of the spin-glass susceptibil- 
ity and concluded that the phase transition does not oc- 
cur. Matsubara et almll reported in the preprint version 
that the same analysis may give a finite transition tem- 
perature. A difference of these two analyses is a tempera- 
ture range of data used in the scaling. We also experience 
this kind of ambiguity in the present paper. 

Determination of good scaling is usually done by the 
looks of scaling plots. All the scaled data are fitted by 
some polynomial functions. One which makes the least 
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standard deviation is selected to be the best scaling re- 
sult. It may change if we use a different set of temper- 
ature/size/time ranges in the analysis. This change is 
systematic. 

For example, we consider a situation when a ratio of ex- 
ponents "f/zv in Eq. I|12|) is larger than the correct value. 
The left-hand-side of the scaling formula becomes smaller 
and smaller as time increases. This tendency is stronger 
at lower temperatures where the correlation time is rel- 
atively large. In order to make scaled data fall onto a 
single curve we must choose t(T) at lower temperatures 
to be smaller than the correct values. This is because 
the universal scaling function \ is a decreasing function 
with respect to t/r(T). Since r(T) oc (T — T sg )~ zu , the 
transition temperature is underestimated when a ratio 
7 / 'zv is larger than the correct value. As we discard 
data at higher temperatures and restrict our analysis to 
lower temperature ranges, the transition temperature is 
more and more underestimated. Contrarily, the transi- 
tion temperature is more and more overestimated when 
the ratio 7/ 1 zv is smaller than the correct value. When 
the ratio is correct, the correct transition temperature is 
obtained and the value is robust against changes of the 
temperature range. 

Our new criterion for the scaling is to find a temper- 
ature and critical exponents so that the obtained results 
are independent from the temperature range of data used 
in the scaling analysis. First, we set one temperature 
range. We obtain the transition temperature which gives 
the best scaling plot at each point of 7 / ' zv. The obtained 
transition temperature is plotted against y/zv. Chang- 
ing the temperature range we perform the same analy- 
sis and obtain another plot for the transition tempera- 
ture. These plots exhibit systematic behavior as men- 
tioned above. They cross at one point, which is the most 
probable estimate for the transition temperature. 

Figure^shows an example of the analysis with the new 
scaling criterion in the ferromagnetic Ising model in two 
dimensions. The lattice size is 499 x 500 with skewed peri- 
odic boundary conditions. The Monte Carlo step is 2400. 
The first 300 steps are discarded, which contain initial re- 
laxation. Simulations are performed at nineteen different 
temperatures ranging 2.276 < T < 2.35. Averages over 
124000 independent Monte Carlo runs are taken. When 
a ratio "//zv is larger, estimates of the transition temper- 
ature go down as the temperature range shrinks to lower 
temperatures. They go up when a ratio is smaller. The 
plots cross near the exact point. Deviation from the exact 
value is larger for the exponent zv as shown in Fig. ^b). 
Crossing behavior is rough in this case. An estimation of 
the exponent is difficult compared to the transition tem- 
perature. We perform this scaling procedure to estimate 
the spin-glass transition temperature and the chiral-glass 
transition temperature. 
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FIG. 1: Results of the finite-time-scaling analysis with our 
new criterion in a ferromagnetic Ising model in two dimen- 
sions. An exact estimate using z = 2.165(10)^ is depicted by 
a symbol. Boxes depict error bars, (a) T c versus "//zv. (b) zv 
versus "//zv. A thin line is an exact line regardless a value of z. 



III. RESULTS 
A. Relaxation of the susceptibility 

First, we observe raw relaxation functions of the glass 
susceptibility. We must check time scale when the finite- 
size effects appear in the relaxation functions. Data must 
be free from the finite-size effects because the noncquilib- 
rium relaxation method is based upon taking the infinite- 
size limit before the equilibrium limit is taken. The sus- 
ceptibility data influenced by the size effects exhibit con- 
verging behavior even though it should diverge in the 
infinite-size limit. It misleads us into thinking that the 
temperature is in the paramagnetic phase. 

Figure [21 shows the appearance of the finite-size ef- 
fects of Xsg(i) and Xc S {t). When a size is small, a 
nonequilibrium relaxation function deviates from the 
size-independent relaxation curve at a characteristic 
time. The system finds itself finite at this time. It be- 
haves as an infinite-size system before the time. There 
is no difference between data of L = 39 and those of 
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FIG. 2: Finite-size effects of relaxation functions of XsgW 
and Xcg(t)- The temperature T = 0.45 is near our es- 
timate for the transition temperature. There is no size 
dependence for L = 39 within 10 Monte Carlo step. 



L = 49. Therefore, relaxation functions of L = 39 within 
t < 10 5 are regarded as of the infinite-size system. Wc 
mainly use this set of size and time limits in this paper. 
Number of random bond configurations is typically more 
than several hundreds. It varies depending on system 
sizes, Monte Carlo steps, and physical quantities. For 
example, quite a few configurations are necessary in or- 
der to obtain beautiful relaxation functions for the Binder 
parameter and the correlation length. 

It is also noted in Fig. [3 that behavior of the finite- 
size effects of the chiral-glass susceptibility is opposite 
to the spin-glass one. The spin-glass susceptibility devi- 
ates from the diverging curve just to converge to a finite 
value. It is natural behavior. On the other hand, the 
chiral-glass susceptibility increases when the size effects 
appear. It eventually converges to a finite value and is 
overtaken by the infinite-size one. This strange behav- 
ior is also observed in the ± J Heisenberg model in three 
dimensions^* 

Figure EJa) shows relaxation functions of Xsg(^) and 
Xcg(t) for T > 0.47. This is a plot on the high temper- 
ature side. It is clear that there is no spin-glass/chiral- 
glass transition at T > 0.50. Both susceptibility exhibit 
converging behavior. They exhibit algebraic divergence 
at T = 0.47. It suggests that the spin-glass transition 
occurs at a temperature near T = 0.47. Figure |3Jb) 
shows a plot at lower temperatures. The spin-glass sus- 
ceptibility diverges algebraically with exponents decreas- 
ing with the temperatures. The exponent is largest at 
T = 0.45. We consider it to be near the transition tem- 
perature. The same behavior is observed in the ±J Ising 
model in three dimensions below the spin-glass transi- 
tion temperaturei 37 ! 38 Slow dynamics appear in the low- 
temperature spin-glass phase. Relaxation of the glass 
susceptibility may become algebraically slow. The chiral- 
glass susceptibility also exhibit algebraic divergence at 
lower temperatures. The amplitude as well as the expo- 
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FIG. 3: Relaxation functions of Xsg and Xcg at various tem- 
peratures, (a) Plots at T > 0.47: high temperature side, (b) 
Plots at T < 0.45: low temperature side, (c) Exponents of 
algebraic divergence in Xsg and Xcg at lower temperatures. 



ncnt is dependent on the temperatures. 

Temperature dependences of the exponents are shown 
in Fig. Etc). Algebraic divergence starts appearing at 
T ~ 0.47. An exponent of the spin-glass susceptibility 
takes a maximum value near T — 0.45. It is consid- 
ered in the vicinity of the transition temperature. The 
exponent decreases as the temperature decreases. This 
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temperature dependence is linear with T: 

j/z(T)v = 0.81T, l/z(T)~0.4T. (13) 

Here, we have used our final estimate for a ratio j/v = 
2.0(2) obtained in Sec. IIII Dl The temperature depen- 
dence of l/z in the Ising model with Gaussian bond dis- 
tributions is reported by Komori et al££ to be 1/ z(T) = 
0.17T. The low-temperature phase in the present XY 
model is similar to that of the Ising model. 

Temperature dependences of an exponent of the chiral- 
glass susceptibility is strange compared with the spin- 
glass one. The algebraic divergence begins at T = 0.47. 
As the temperature decreases, the exponent first in- 
creases and takes a maximum value at T = 0.42 and 
then decreases with temperature. The temperature of 
the maximum exponent is lower than that of the spin- 
glass one. 

From the relaxation functions of the susceptibility it is 
very possible that the spin-glass phase transition occurs 
at T = 0.45 ± 0.03. The chiral-glass transition also oc- 
curs within this temperature range. Since the chirality 
trivially freezes if the spin freezes, these two transitions 
may occur simultaneously. The low-temperature relax- 
ation functions of the spin-glass susceptibility exhibit a 
property of the spin-glass phase: an inverse of the dy- 
namic exponent is linearly dependent on the tempera- 
ture. We check these findings performing the finite-time 
scaling analysis. Stress is put upon consistency of a tran- 
sition temperature and a ratio of critical exponents ^ j zv 
at that temperature. 



B. Finite-time-scaling analyses 

We perform a finite-time-scaling analysis with our new 
criterion. Figures Ufa) and^Jb) show results of the scal- 
ing on the spin-glass transition. For each value of ^ / zv 
we choose T sg and zv so that fitting error of the scaled 
data becomes smallest. The fitting is done by a polyno- 
mial expression up to the 7th order. Relaxation data of 
the first two thousand steps are discarded. We considered 
them to be the initial relaxation. The scaling results are 
dependent on the temperature range. The dependence is 
also as described in Sec. Ill Al when 7 / zv is larger, the 
transition temperature becomes underestimated as high 
temperature data are discarded, and vice versa. From 
the crossing point we obtain our estimates as 

T sg = 0.455 ±0.015 (14) 
j/zu = 0.355 ± 0.005 (15) 
zv = 4.8 ±0.4. (16) 

Figure 0Jc) shows a scaling plot using the obtained re- 
sults. All data fall onto a single scaling function very 
nicely. The obtained transition temperature is higher 
than other estimates. Granato^i gave T sg = 0.39(2), 
v = 1.2(2) and z = 4.4(3). A value of zv is consistent 



with our value. Lee and YoungSl gave for the Gaussian 
model T sg = 0.34(2) and u = 1.2(2). Our estimate for 
v is v = 0.85(15), where we use a value of z which will 
be obtained in Sec. IrTlDl fEa. QQ) ). Though the transi- 
tion temperature is quite different, critical exponents are 
consistent. 

Values of T sg and 7/ 1 zv obtained by the scaling anal- 
ysis agree with the relaxation function of the suscepti- 
bility in the previous subsection. As shown in Fig. |2f c ) 
an exponent of the spin-glass susceptibility takes a max- 
imum value at around T — 0.45 and the exponent is 
j/zv = 0.348. Both values are consistent with the scal- 
ing results. A scaling estimation of T sg = 0.4 is possible if 
we artificially choose 7/ 'zv and the temperature range of 
the data used in the scaling analysis. An example is a set 
of choices: 7 / zv = 0.366 and T = 0.50 ~ 0.56 as shown 
in Fig. 0Ja) with circles. However, it is inconsistent with 
the relaxation functions of the spin-glass susceptibility. 
It does not diverge with an exponent j/zv = 0.366 at 
T = 0.4. Only the scaling results that are independent 
from the temperature ranges (a crossing point in the fig- 
ure) become consistent with the raw relaxation functions. 

Figures |3a),EIb) and^c) shows results on the chiral- 
glass transition. Estimations of the transition tempera- 
ture and critical exponents are 

T cg = 0.467 ±0.010 (17) 
J K /z K v K = 0.30 ±0.01 (18) 
z K v K = 4.7 ±0.2. (19) 

Subscripts k of the critical exponents denote chirality. 
Kawamura and Li gave T cg = 0.39(3), v K = 1.2(2), 
7k/^k = 1.85(20), and z K = 7.4(10). Our estimation for 
the transition temperature is higher than their value. We 
will obtain a value of z K = 6.3(5) from a relaxation func- 
tion of the chiral-glass correlation length in Sec. 1111 Dl 
Then, an exponent v K = 0.75(10) is obtained. This value 
is also smaller than their value. However, a ratio of ex- 
ponents, j K /v K = 1.9(2), is consistent. As observed in 
Fig. IsJfc) an exponent j K /z K v K of the chiral-glass suscep- 
tibility takes a maximum value of ^ K /z K v K = 0.315 at 
T = 0.42. This temperature is consistent with the esti- 
mation of Kawamura and Li. However, this set of T and 
7 K / ' z k v k is inconsistent with the scaling results of Fig. 
|3Ja). To the contrary, our scaling estimations are con- 
sistent with the relaxation functions of the chiral-glass 
susceptibility. We consider that the chiral-glass transi- 
tion occurs at T ~ 0.47. Algebraic divergence begins at 
this temperature. The temperature of the exponent max- 
imum does not correspond to the transition temperature 
of the chirality. 

The spin-glass transition temperature and the chiral- 
glass transition temperature agree well within the error 
bars. Therefore, both transitions are considered to occur 
simultaneously Estimations are consistent with the raw 
relaxation functions. The consistency guarantees that 
our analyses are performed correctly. 
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FIG. 4: Results of the present finite-time-scaling analysis 
on the spin-glass transition, (a) T sg versus y/zv. (b) zv 
versus y/zv. (c) A scaling plot using the obtained parameters. 



FIG. 5: Results of the present finite-time-scaling analysis 
on the chiral-glass transition, (a) T sg versus y/zv. (b) zv 
versus y/zv. (c) A scaling plot using the obtained parameters. 



C. Relaxation of the Binder parameter 

In order for another check of our scaling results and 
for obtaining the dynamic exponent z independently we 
have calculated the noncquilibrium relaxation functions 
of the Binder parameter at the estimated transition tem- 
perature. 

Figure EJa) shows relaxation functions of the Binder 
parameter of the spin, g ss (t), for various system sizes. 



When a size is small, a relaxation function deviate from 
the sizc-indcpcndcnt relaxation curve. Appearance of the 
finite-size effect is normal as in the case of the spin-glass 
susceptibility It just converges to a finite value when 
the size effect appears. The size-independent curve ex- 
hibits an algebraic divergence with an exponent d/z with 
z ~ 6.0 ±0.5. We have also obtained the relaxation func- 
tions at T = 0.42 and T = 0.47. (Figures not shown in 
this paper.) These temperatures are within the error bar 
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FIG. 6: Relaxation functions of the Binder parame- 
ters at T = 0.45. (a) A log-log plot for the spin- 
glass one. Algebraic divergence is observed, (b) A semi- 
log plot for the chiral-glass one. No divergence is ob- 
served. When the size effects appear, the Binder parame- 
ters take negative values first and converge to positive values. 



of the spin-glass transition temperature. The algebraic 
divergence of the Binder parameter is also observed. It is 
noted that the dynamic exponent z is systematically de- 
pendent on the temperature. The values are z = 6.5±0.7 
at T = 0.42 and z = 5.7 ± 0.5 at T = 0.47. They are 
consistent with the temperature dependence of z in Eq. 
(|13[l : 1/ z(T) ~ 0.4T within error bars. The dynamic ex- 
ponent is estimated to be z = 6.0 ± 0.8 within an error 
bar of the spin-glass transition temperature. 

Relaxation functions of the Binder parameter of the 
chiral-glass shows strange behavior. The algebraic diver- 
gence as in the spin-glass one cannot be observed at all in 
Fig. EJb) . A value of the Binder parameter remains zero 
before the finite-size effects appear. Then, it decreases 
and takes negative values. Taking a minimum value it 
increases with time and eventually converges to a posi- 
tive value, which is consistent with a result of Kawamura 
and Li^i. This strange finite-size effect possibly has the 
same origin with that of the chiral-glass susceptibility as 
shown in Fig. [5] The chiral-glass Binder parameter is 
considered to take zero value at T > T cs in the infinite- 
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FIG. 7: Estimations of a correlation length £ Bg (t) by fitting 
a correlation function f B g(l). Data at t = 10 2 are depicted 
by circles and those at t = 10 5 are depicted by squares. 



size system. 



D. Relaxation of the correlation length 

In this subsection the final check is made on our anal- 
ysis concerning the spin-glass transition. We also answer 
questions and claims regarding use of the nonequilibrium 
relaxation method in spin-glass problems. 

Figure \7\ shows the spin-glass correlation functions at 
Monte Carlo steps £ = 10 2 and t = 10 5 . The temperature 
is near the spin-glass transition temperature. We have 
taken averages over three directions, x, y and z, between 
spins. The distance is limited within L/2 because we 
have imposed skewed periodic boundary conditions. Our 
simulations start from the perfectly random spin config- 
uration. There is no correlation between spins and the 
correlation length is zero at t = 0. The spin-glass cor- 
relation grows with time. The correlation decays fast at 
t = 10 2 but decay is slow at t = 10 5 . We estimate the 
correlation length fitting the correlation functions by an 
expression Eq. lllj l. Here, we have set 7/1/ = 2.4. We 
obtain a ratio of j/v ~ 2.1(3) from the scaling result 
of 7/ ' zv = 0.355(5) together with an estimate of the dy- 
namic exponent z = 6.0(8) by the Binder parameter. We 
performed the fitting with several values of 7/1/ in this 
range. Amplitudes of the correlation lengths change with 
7/V. When a value of 7/z/ is small, £ sg is overestimated: 
it becomes 40% larger if we set 7/1/ = 1.8. However, the 
dynamic-scaling behavior, £ ~ t 1 / 2 , and the exponent z 
remain the same. We discard data of the correlation func- 
tions which are lower than 10 -3 in the fitting procedure. 
This is resolution of the present simulations. Numerical 
data may take negative values below this resolution limit 
even though the definition Eq. Q demands positive val- 
ues. The obtained correlation length is not influenced by 
such garbage data. The spin-glass correlation length at 
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FIG. 8: Relaxation functions of (a) the spin-glass cor- 
relation length and (b) the chiral-glass correlation length. 
Linear sizes of lattices are denoted by L. Fitting lines 
for the size-independent data are £ sg (f) = 0.4t 1//5,1 and 
?c g (i) = 0.3t 1//6 ' 3 , respectively. The size effect appears 
when the correlation length reaches L/10 in every case. 



t = 10 2 is £ sg = 1.00, and £ sg = 3.84 at t = 10 5 . The 
growth is very slow. We have also obtained the chiral- 
glass correlation length in the same manner. 

Figure [H] shows relaxation functions of the spin-glass 
correlation length and the chiral-glass correlation length. 
At each Monte Carlo step we calculate the correlation 
functions and obtain the correlation lengths. The cor- 
relation length exhibits algebraically diverging behavior 
until t = 10 6 . There is no size dependence between data 
of L = 39 and those of L = 49. We can regard them the 
infinite-size data. Random configuration number of the 
L = 39 data is 235, and that of the L = 49 data is 122. 
Deviations from the fitting line arc larger for the L = 49 
data. The dynamic exponent of the spin-glass correla- 
tion length is2 = 5.1±0.5. This value is a little smaller 
than the estimation by the Binder parameter, which gave 
z = 6.0 ±0.8. Both estimates are independent. Thus, we 
take an average of them and obtain our final estimate. 



z = 5.6 ± 0.6 (spin) 
z K = 6.3 ± 0.5 (chirality) 



(20) 
(21) 



Using the value of spin and j/zis = 0.355(5) we obtain a 
ratio of the critical exponents of the spin-glass suscepti- 
bility as 7/1/ = 2.0(2). 

It is remarkable that the spin-glass correlation length 
starts the diverging behavior almost from the first step. 
A correlation length is about 0.5 lattice spacings. The 
chiral-glass correlation length also starts the diverging 
behavior at t = 50 where a correlation length is also 0.5 
lattice spacings. This is a characteristic length above 
which we can observe the algebraic divergence of the cor- 
relation length. It is considered to be determined by 
resolution of simulations. The figure also shows how the 
finite-size effects appear in the relaxation functions of 
the correlation lengths. In both plots the finite-size ef- 
fects appear when the correlation lengths reach L/10. 
We have observed the size effect at other temperatures. 
It also appears when £ ~ L/10. 

The size effects enhance the correlation length. The 
skewed periodic boundary conditions are imposed in this 
paper. Boundary effects produce correlation of a spin 
with itself. Therefore, the correlation is always overesti- 
mated by the size effects. Then, the correlation length is 
also overestimated. The spin-glass correlation length of 
L = 9 reaches £ sg = 6 in the equilibrium limit as shown 
in Fig. IHIa). In a lattice of L = 9 the largest distance 
is L/2 = 4.5. The correlation length exceeds this limit. 
Therefore, such a value does not have a physical mean- 
ing at all. This is a consequence of the size effect. One 
may think that the size effect appears when the correla- 
tion length reaches the system size. However, it appears 
when £ reaches L/10. Finite-size effects are very strong 
in spin glasses. We must pay much attention to the size 
effects on numerical data. 

These lines of evidence may answer the question why 
the nonequilibrium relaxation method can handle the 
spin-glass problems. Relaxation processes in the scaling 
region know the critical phenomenon from the very early 
time steps. There is a characteristic length that we can 
observe the critical behavior. This is £ = 0.5 lattice spac- 
ings in the present simulations. This length is surpris- 
ingly small. This is because resolution of our simulations 
is very sharp. It is an advantage of the nonequilibrium 
relaxation method. Critical relaxation continues until 
the finite-size effects appear. It is a time when £ reaches 
L/10. If the temperature is off the transition temper- 
ature, critical relaxation continues until the correlation 
length reaches its equilibrium finite value at the temper- 
ature. What is slow in spin glasses is a relaxation process 
after the finite-size effects appear. For example, £ sg of 
L = 9 in Fig. 02a) deviate from the size-independent line 
at t ~ 30. It reaches the equilibrium limit at t = 10 5 . 
Computations before this time are discarded in the equi- 
librium simulations. On the other hand, those before the 
size effects appear (t ~ 30) are utilized in the nonequi- 
librium relaxation method. Therefore, the critical phe- 
nomenon of spin glasses is easily observed if we take the 
nonequilibrium relaxation approach. Another advantage 
of this approach is that the system size can be considered 
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infinity. The size effects arc very strong and strange in 
spin glasses as shown in this paper. We can get rid of 
these effects by this method. These are the reasons why 
the nonequilibrium relaxation method is successful in the 
spin-glass problems. 



IV. SUMMARY AND DISCUSSION 

In this paper we made clear that the simultaneous spin- 
glass and chiral-glass transition occurs in the ±J XY 
model in three dimensions. We have eliminated ambi- 
guity of the scaling analysis by introducing a new crite- 
rion: correct results arc independent of the temperature 
range of data used in the scaling. The results are checked 
and confirmed by the relaxation functions of the spin- 
glass susceptibility, the Binder parameter and the spin- 
glass correlation length. They are consistent with each 
other. The evaluated transition temperature is higher 
than other estimationsj2L2i while the critical exponents 
are consistent. 

The transition temperature of the spin-glass transition 
and that of the chiral-glass transition agree well within 
the error bars. A ratio of the static critical exponents, 
7/1^, takes almost the same value between two. They are 
•y/i/ = 2.0(2) for the spin and 7 K /^ K = 1-9(2) for the chi- 
rality. This agreement suggests that two transitions are 
qualitatively equivalent. The dynamic exponent z of the 
spin-glass obtained by the Binder parameter and that ob- 
tained by the correlation length becomes consistent when 
T sg = 0.47. This temperature is an upper bound of the 
error bar and is very close to the mean value of the chiral- 
glass transition temperature. 

The spin-glass susceptibility and the Binder parameter 
exhibit normal behavior of the second-order phase transi- 
tion. Both exhibit algebraic divergence at the transition 
temperature. An exponent of the divergence is strongest 
at the transition temperature. Contrarily, those of the 
chiral-glass transition show strange behavior. The chiral- 
glass susceptibility diverges algebraically at the transi- 
tion temperature. However, an exponent of the suscep- 
tibility takes a maximum value at a lower temperature. 
Relaxation of the Binder parameter shows no divergence. 
It just remains zero in the infinite-size system. If the 
transition is driven by the chirality degrees of freedom, 
the Binder parameter should exhibit a kind of critical 
behavior as the spin-glass one exhibits. Relaxation func- 
tions of the chiral-glass correlation length exhibit the 
same behavior as the spin-glass one. They begin the 
algebraic divergence when £ S g/cg ~ 0-5 lattice spacings 
and stop it at £ S g/cg ~ L/10 when the finite-size effects 
appear. The size effects appear later in the chiral-glass 
because £ cg < £ sg . Kawamura and Li2i commented that 
the spin-chiral separation occurs at T = T cg at t = 10 s 
Monte Carlo steps in the standard Metropolis dynam- 
ics. We have carried out simulations up to t = 10 6 steps 
and found no sign of the separation in behaviors of the 
infinite-size system. 



The lines of evidence can be explained if one considers 
that the glass transition is driven by the spin degrees of 
freedom. The chirality is defined by spin variables. It 
trivially freezes if the spin freezes. A critical property 
may be observed by some quantities of the chirality, but 
it may not in other quantities. The former example is the 
chiral-glass susceptibility and the latter one is the Binder 
parameter. This is because the chirality is the secondary 
degrees of freedom in this phase transition. Every analy- 
sis and the results are consistent in the spin-glass transi- 
tion. This is a strong support for our argument that the 
transition is driven by the spin. 

The present analyses are based upon numerical simu- 
lations. The results always have error bars. We cannot 
prove that the spin-glass transition temperature is ex- 
actly same with the chiral-glass one. Our estimations 
for the transition temperatures agree within error bars. 
However, the mean value of the chiral-glass transition 
temperature is a little larger. In order to make a dis- 
tinction between two transition temperatures we must 
carry out simulations 400000 times as large as the present 
ones. The estimation is as follows. If the mean val- 
ues are correct, Xsg(t) should exhibit converging behav- 
ior when T = 0.467, while Xcg(^) remains diverging. 
We observed converging behavior of Xsg(i) at t = 10 4 
when T = 0.56 as shown in Fig. |3{a). Therefore, the 
converging behavior is observed at t = 3.3 • 10 8 when 
T = 0.467 because r cx (T — 0.455)" 4 8 . The size ef- 
fects appear when £ ~ L/10 and £ sg (i) = OAt 1 / 5 - 6 . 
Then, a linear lattice size L = 130 is necessary to ob- 
tain size-independent data until t = 3.3 • 10 s . Therefore, 
3.3- 10 3 x (130/39) 3 ~ 400000 times larger simulations are 
necessary. It is almost impossible to make the distinction 
by raw relaxation functions. 

Why the spin-glass transition has not been observed 
until quite recently? The first reasonable answer we con- 
sider is strong finite-size effects. As shown in Fig. [S] the 
size effects appear when the correlation length reaches 
L/10. Physically-meaningful spin-glass correlation func- 
tions are restricted within this length scale. Those out- 
side this scale are plagued by the size effects. Since the 
susceptibility collects all the correlation functions, it is 
also under the strong influence of the size effects. Most 
of information extracted from the susceptibility may be 
unphysical one if the lattice size is small. Observations 
of the spin-glass susceptibility by equilibrium simulations 
with small lattices may suffer from this difficulty. The 
second reason is the boundary conditions. Most simu- 
lations employ the periodic boundary conditions. They 
may produce a domain-wall spin state, or a spiral spin 
state in continuous spin models. Particularly, when the 
system possesses frustration, nontrivial states may be 
produced. On the other hand, the nonequilibrium re- 
laxation method achieves the infinite-size system. The 
obtained data are always independent of the finite-size 
effects and the boundary conditions. True physical prop- 
erties are easily extracted. 

Recently, the spin-glass transitions are observed by a 
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finite-size scaling analysis on the spin-glass correlation 
length^ The quantity is considered to be less influenced 
by the size effects. In the present paper we have also ob- 
tained relaxation functions of the correlation length. It 
is evaluated directly from the correlation functions. Un- 
physical garbage data of correlations are not contained 
in our estimations. Therefore, the spin-glass transition 
is clearly exhibited. It is also found that the correla- 
tion length exhibits critical behavior from the first step. 
Analyses on the correlation length will become a stan- 
dard approach in the spin glass problems. 
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